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Abstract 



In this paper, we present a realistic application of the coherence pro- 
tection method proposed in the previous article. A qubit of information 
encoded on the two spin states of a Rubidium isotope is protected from 
the action of electric and magnetic fields. 

1 Introduction 

The coherence protection method presented in the previous article is applied to 
a Rubidium isotope, in which the information part corresponds to the spin of 
the exterior electron, whereas the orbital part of the wavefunction plays the role 
of the ancilla. 

The different steps of our protection method are illustrated on this specific 
example. Adding the ancilla is achieved through pumping the atom from the 
level 5s into the shell 60f. The coding and decoding matrices are applied through 
the non-holonomic control technique. Finally, the protection step is achieved 
through the simultaneous emissions of three photons : two emissions are stimu- 
lated, whereas the third is spontaneous and constitutes the irreversible process 
needed in this step. 

The paper is organized as follows. In the second section, we present the 
system and motivate our choice. In the third section we review each step of 
our method in detail and show their implementations on this specific physical 
configuration. Moreover, experimental problems and limitations are discussed. 

2 Presentation of the system 

In this paper, we show how to protect one qubit of information encoded on the 
two spin states of the ground level 5s of the radioactive isotope 78 Rb against the 
action of M — 6 error- inducing Hamiltonians E m . For numerical calculations 
we considered 3 magnetic Hamiltonians 



Let us first motivate the choice of the Rubidium atom. Alkali atoms like Rb 
are very interesting for our purpose because of their hydrogen-like behavior : 
indeed, such an atom is the "natural" compound of an information subsystem, 
the spin part of the wavefunction, and an ancilla, the orbital part of the quantum 
state. As we shall see, one can easily increase the dimensionality of the ancilla by 
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Figure 1: Spectrum of 78 Rb: The useful part of the spectrum of Rubidium is 
represented. 



simply pumping the atom towards a shell of higher orbital angular momentum 
L. 

Among all alkali systems we chose 78 Rb because of its spectroscopic features 
(FigCJ: in particular, 78 Rb has no hyperfine structure (its nuclear spin is 0) 
which ensures that the ground level 5s is degenerate (this is necessary for the 
projection scheme as we shall see below). Moreover it has a long enough lifetime 
(r ~ 17.66 min) for the proposed experiment. 

3 The different steps of the method 

We shall now review each step of our method in detail. The information we want 
to protect is initially encoded on the two spin states = \5s,j — ^,rrij = — 5) 
and ^2) = 5s, j — ^,nij = \) of the ground level 5s of the atom: the two- 
dimensional (/ = 2) information space "Hi — Span , \v2)\ is spanned by 
these two states. The first step of our scheme consists in adding an ancilla A to 
the information system. In the present setting, the role of A is played by the 
orbital part of the wavefunction. In the ground state (L = 0), its dimension is 
A = 2L + 1 = 1 (roughly speaking, there is no ancilla). If we want to protect 
one qubit of information against M = 6 error-inducing Hamiltonians, we have 
to increase the dimensionality of the ancilla up to A = M + 1 = 7 (according 
to the Hamming bound presented in the previous paper), by pumping the atom 
towards a shell nf (L = 3). We choose the highly excited Rydberg state 60/ 
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so as to make the fine structure as weak as possible (the splitting for 60/ is 
approximately 10 _5 cm _1 ), which shall be neglected in a first approach : thus, 
heN = IxA = 2x7 = 14 basis vectors of the total Hilbert space Tt = Hi^Ti-A 
are almost perfectly degenerate ; the validity of this approximation will be 
discussed at the end of this section. To be more specific, the pumping is done 
in the following way 
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In other words, using the terminology of the previous sections, the information 
initially stored in TLi is transferred into 



C = Span 
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The choice of the subspace C might appear arbitrary at this stage, but it will be 
justified later by the practical feasibility of the projection process onto C. Note 
that C is an "entangled" subspace, whose basis vectors {|7i)} i=1 2 are generic 
entangled states of the spin and orbital parts: this means that the projection 
step will not consist in a simple measurement of the ancilla but will involve a 
more intricate process we shall describe in detail later. 

Practically, the pumping can be achieved as follows. Three lasers are ap- 
plied to the atom: the first laser is right polarized and slightly detuned from the 
transition (5s «-> 5p) whereas the second and third lasers are left polarized and 

slightly detuned from the transitions (^>p3 <-> 5<^|) and ^5d| <-> 60/^j respec- 
tively. The role of the dctunings is to forbid real one-photon processes: thus, 
the atom can only absorb three photons simultaneously and is thereby excited 
from the ground level 5s to the Rydberg level 60/. By using selection rules, one 
can construct the allowed paths represented on Fig|2 these paths only couple 
and 1 1/2) to 1 71) and I72), respectively. 



The second step consists in imposing the coding matrix to the system by the 
non-holonomic control technique: to this end, we submit the atom to nc = 34 
control pulses of timings {U} i=1 34 , during which two different combinations of 
magnetic and Raman electric Hamiltonians are alternately applied (see Fig|3J). 
To be more explicit, during odd- numbered pulses ("A" type pulses) we apply 
the constant magnetic field 

_^ / B X = 7 10- 3 T \ 
B = B y = 8.2 10- 3 T 

\ B z = -6.8 10~ 3 T / 

associated with the Zeeman Hamiltonian Wz, and two sinusoidal electric laser 
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Figure 2: Ancilla adding by Pumping. Photon polarization and involved sub- 
Zeeman levels are represented. 
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whose frequencies lor and u>' R are respectively slightly detuned from the two 
transitions (60/ 5d,j = |) and (60/ 5d,j = f ) (detunings <5 and 5'). The 
characteristic values of these fields are 
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The intensity of the laser beams are typically of the order of 2 10 8 W.cm 1 . The 
Raman Hamiltonian associated with these fields is denoted by Wr^a and the 
total perturbation is P a = Wz + Wr.a- During even- numbered pulses ("B" 
type pulses), we apply the same magnetic field as for A type pulses, which is 
experimentally convenient, and two sinusoidal electric laser fields 
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whose frequencies are the same as above and whose characteristic values are 

= -5.2 lO^.m" 1 
= 8.5 lO'V.m" 1 
= 2.3. 

The Raman Hamiltonian associated with these fields is denoted by Wr^b- The 
corresponding perturbation is P b = Wz+Wr^b- As the fine structure of the level 
60/ is neglected, the unperturbed Hamiltonian Hq is and the total Hamiltonian 
has the following form: H a = P a during "A" pulses, H b = P b during "B" pulses. 
The 34 different timings have been calculated so that 



U{ti, ■■■,t 34 ) = e- lHBr ^e' lHAT33 . . . e~ iHATl 
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Figure 3: Coding step through the non-holonomic control technique. The two 
Hamiltonians H a and Hb are alternately applied to the system during pulses 
of timings {n(ns)} ={3.9763, 6.4748, 4.2274, 3.6259, 2.8717, 3.6281, 7.2263, 
6.4260, 4.8070, 5.0394, 6.5242, 4.8890, 4.2400, 7.3834, 4.8653, 5.4799, 4.5341, 
4.3099, 6.2959, 3.7346, 6.5293, 6.8586, 6.0749, 5.1213, 4.6806, 3.4985, 3.9909, 
4.6701, 4.5168, 6.4702, 4.7787, 5.3476, 3.4567, 3.8009}. The frequencies of the 
laser fields involved in the encoding step are represented on the spectrum of the 
Rubidium atom. The fine structure of the Rydberg level 60/ is not represented. 



checks the correction conditions presented in the previous paper. At the end 
of the coding step, the information is transferred into the code space C = CC, 

encoded on the codewords I hi) = C hi) > 

I J i=l,2 

As one can see on Fig|3| the total duration of a control period (~ 125ns) is 
approximately 10 3 times shorter than the lifetime of 60/ Rydberg state which 
is approximately 0.115rns, and the different pulse timings range between 2.9ns 
and 7.4ns, which are feasible. 

After a short time, due to the action of the error Hamiltonians, the infor- 
mation stored in the system acquires a small erroneous component, which is 
orthogonal to the code space C. Then, we decode the information through the 
application of the matrix C -1 . To this end, we reverse B and the detunings S 
and 5' ', while leaving all the other values unchanged (this amounts to taking the 
opposite of Hamiltonians H a and and apply the same sequence of control 
pulses backwards: we start with an "A" pulse whose timing is T34, then apply 
a "B" pulse during t 33 , etc. (see FigQJ. The decoding step yields an erroneous 
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Figure 4: Decoding step through the non-holonomic control technique. We 
reverse the magnetic field and the detunings of electric fields, as represented on 
the spectrum of the Rubidium atom, and apply the same control sequence as 
for coding (same timings) in the reverse way. 



state whose projection onto C is the initial information state. 

In the last step, we project the erroneous state vector onto the subspace C 
to recover the initial information. Projection is a non-unitary process which 
cannot be achieved through a Hamiltonian process, but requires the intro- 
duction of irreversibility. To this end, we make use of a path which is sym- 
metric with the pumping step, consisting in two stimulated and one sponta- 
neous emissions. To be more explicit, we apply two left circularly polarized 
lasers (see FigEJ) slightly detuned from the transitions (60/ < — ► 5d,j — §) and 
(5d, j — | < — > 5p, j = |j. Due to these laser fields, the atom is likely to fall to- 
wards the ground state and emit two stimulated and one spontaneous photons. 

If a circularly right-polarized spontaneous photon is emitted, the selection 
rules show that the only states to be coupled to the ground level are I71) and 
I72) to \v\) and |z^), respectively (see Fig|3J). In other terms, the emission of 
such a spontaneous photon brings the "correct" part of the state vector back 
into Hi = Span , \v2)\. On the contrary, the other cases - "left-polarized", 
"linearly-polarized spontaneous photon", or "no photon at all"- do not lead to 
the right projection process. 

The "left-polarized photon" and "no photon emitted" cases are quite un- 
likely, since their probability is proportional to the square of the error ampli- 
tude, that is to the square of the very short Zeno interval. By contrast, the 
"linearly polarized photon" case is quite annoying because it mixes the two 
paths 1 71) — ► and I72) — > ^2) : one has to get rid of this parasitic pro- 
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Figure 5: Projection path. The lasers involved are marked by solid arrows, the 
spontaneous photon is represented by a dashed arrow. The different polariza- 
tions are specified. The fine structure of the level 60/ is not represented. 
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cess by minimizing its relative probability with respect to the process followed 
by the "right-polarized" photon emission. This can be achieved by launching the 
78 Rb atom, previously cooled, into a Fabry-Perot cavity, in an atomic fountain 
manner (fine tuning of the lasers driving the 60/ — 5c? and bd — bp transition will 
be necessary to avoid reflection of the external laser radiation from the cavity). 
The decay rate for the 3-photon transition \"fi) — ► \vi) is 



- liVi 



2tt 



fiAi 



d\^ k E 2 



h(A 1 + A 2 ) 



" VkVi e R 



where k s denotes the wave vector of the spontaneously emitted photon, e r 
is the left-polarized photon polarization unit vector, g ^ k s ^ is the density of 

states (normalized to the cavity volume) for the cavity field at k s , and the bar 
denotes averaging over the directions of k s . The transition dipole moments are 
denoted by d a b: during the projective process the states coupled to I71) and I72) 
are 
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respectively. 

The presence of the cavity enhances the density of states for the modes 
propagating paraxially to the z-axis and ensures that 
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where 7 is the decay rate of \5p,j = §,?rij = +|) into \5s,j = \,mj = +5), 
so that the undesired process followed by the 7r-photon emission is relatively 
less important than it were in free space. The dynamics of the density matrix 
elements p a b is governed by the following system (i = 0, 1): 



Pan ~ 



'-■yiVi P~li~1i 1 
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To avoid dephasing which would corrupt the information, the coherence matrix 
clement /0 7l72 must be transferred with the maximum efficiency r\ 
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into p UlU2 . According to the Wigner-Eckart theorem, we have 



- Tit/1 

- 72^2 



^5/2 -3/2 
°3/2 -1/2 1 



^3/2 -1/2 
-l°3/2 1/2 1 



^3/2 1/2 
-l°l/2 -1/2 1 1 



r 5/2 -1/2 
°3/2 1/2 1 



^,3/2 1/2 ^,3/2 3/2 

-l°3/2 3/2 1 -l°l/2 1/2 1 1 



whence 77 = 12\/2/17 « 0.99827. In other words, the probability of error 
during the Zeno projection stage due to the small difference of the Clebsch- 
Gordan coefficient products for the two paths is equal to or less than 1 — 77 w 
0.00173 (the equality is reached if the initial state is (|0) ± \l))/y/2). Note 
that the states 60/, 5d, and 5p have finite lifetimes (see Fig^l. Thus the 



transition rates T liUi must be much larger than 1/rgo/, 



/T 5d , and 



riAi 



/r^p, in order to minimize errors caused by the decay of 



h(A 1 +A 2 ) 

these unstable states. 

To complete the projection step, one has now to transfer the atom back 
into its coherent superposition on the shell 60/: this is achieved by the same 
pumping sequence as in the first step. The mismatch of the Clebsch-Gordan 
coefficient products will cause again the error probability 1 — 77. The information 
is then restored with very high probability and the system is ready to undergo 
a new protection cycle. 



Until now, we have neglected the fine structure splitting of the level 60/, 
which is approximately 2.10~ 5 cm~ 1 and corresponds to a period tj ~ 1.5/zs. 
To conclude this section, we shall now take it into account and see how it affects 
each step of our scheme. 

Obviously the pumping and projection steps will not be affected by the fine 
structure, since the information-carrying vectors {I71) , I72)} belong to the same 
multiplet (J = 5/2). 

The coding and decoding steps are neither modified by the existence of 
the fine structure. Since the typical period of the fine structure Hamiltonian 
tj ~ 1.5 [is is more than 10 times longer than the total duration of the coding 
or decoding steps, it is legitimate to neglect its effect. 

The influence of the fine structure on the free evolution period during which 
errors are likely to occur is more complicated to study in the general case. 
However, two simple limiting regimes can be considered. If the spectrum of the 
coupling functions fm(t)'s is very narrow (i.e. if the variation timescale of the 
fm(t)'s is much longer than 17), one can show that our scheme applies directly 

as though there were no fine structure, provided the error Hamiltonians |-E m | 

are replaced by , where is obtained from E m by simply setting to 

zero the rectangular submatrices which couple the two multiplets (J = 5/2, 7/2). 
The second limiting regime corresponds to a very broad spectrum for the /m(i)'s 
(variation timescale much shorter than Tf): in that case, one can show that our 
scheme applies provided one chooses a Zeno interval multiple of 77. 
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Finally, it must be emphasized that Rydberg atoms, though long-lived are 
also very sensitive to collisional processes as well as Doppler effects, which result 
in a fast coherence loss ; nevertheless, we hope that in a single (or few) atom 
experiment using cold atoms the coherence lifetime can be extended to hundreds 
of ns. These effects have been omitted in all this section in which we just 
intended to provide a pedagogical demonstration of our method on a simple 
physical system ; but an experimentally feasible setup should obviously deal 
with these unavoidable drawbacks of Rydberg states. 

4 Conclusion 

In this paper, a realistic application of our coherence protection method has been 
proposed : it has been shown that, in principle, one qubit of information encoded 
on the spin states of a Rubidium isotope can be protected from the action 
of parasitic electric and magnetic fields. The different steps of our technique 
can be implemented on this specific example : adding the ancilla is achieved 
through pumping ; information is coded and decoded through non-holonomic 
control ; projection is achieved by a three-photon process, involving spontaneous 
emission. 

Practical feasibility of our scheme has been discussed, and experimental 
problems have been raised. However relevant, these limitations do not restrict 
the applicability of our method, and the pedagogical example considered here 
demonstrates its implementability. 
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